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A dynamical pair breaking effect is evidenced at very low excitation energies. For this purpose, a new set of 
time-dependent coupled channel equations for pair-breaking in superfluid systems are deduced from the variational 
principle. These equations give the probability to destroy or to create a Cooper pair under the action of some pertur- 
bations or when the mean field varies in time. The odd-even effect in fission is investigated within the model as an 
example. For this purpose, the time-dependent probability to find the system in a seniority-one or in a seniority-two 
state is restricted in the sense that the perturbations are considered only in the avoided crossing regions. 

Key words: Pair breaking effect. Fission, Landau-Zener effect 
PACS: 21.10.Pc, 24.10.Eq, 24.10.-i, 25.85. -w 



1. Introduction 

In atomic or nuclear molecules, fermions move in a 
common field generated by several centers of force. The 
potential depends on the so-called generalized coordi- 
nates. The most used generalized variables are function 
of the inter-nuclear distances and their angular depen- 
dencies. In such a molecular description, the motion 
of the nuclear centers is considered adiabatically slow 
compared with the rearrangement of the mean field. The 
atoms or the nuclei share their valence electrons and 
their outer-bound nucleons, respectively. Enhancements 
and structures in the energy dependence of the cross sec- 
tion m atomic and nuclear [3, 4] reactions reflect a 
special signature of the presence of molecular orbitals. 
The behavior of the cross section is usually explained in 
terms of the Landau-Zener effect. 

The Landau-Zener effect 0] describes the non- 
adiabatic transitions at avoided crossing regions be- 
tween potential curves f?, '81 or energy levels 
Two levels with the same good quantum numbers as- 
sociated to some symmetries cannot intersect and ex- 
hibit avoided level crossing regions. A small pertur- 
bation energy is always available in these regions, be- 
ing responsible for the level slippage mechanism. Re- 
cently, the time-dependent equations that describe quan- 
titatively the Landau-Zener promotion mechanism were 
generalized for molecular potentials that include pairing 
residual interactions 1 11]. As evidenced below, this gen- 
eralization also provides a way to formulate a theory for 
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a new dynamical pair breaking mechanism. Originally, 
the pairing formalism, reflected in the BCS equations, 
was used to explain superconducting states in terms of 
two electrons pairs with opposite momenta and spin 
near the Fermi surface ll2ll . 



2. Formalism 

The dynamical pair breaking effect emerges from 
a new set of coupled channel equations deduced for 
the time-dependent probability to find the system in a 
seniority-one state or in a seniority-two one. The varia- 
tion of the mean field is considered slow enough, so that 
fermions follow eigenstates of the instantaneous mean 
field. In such an approximation, if the interactions pro- 
duced in the avoided crossing regions or those due to the 
Coriolis coupling are not taken into consideration, other 
perturbations between two different states are not possi- 
ble. In the following, the calculations are restricted only 
for perturbations produced in the avoided crossing re- 
gions. This approximation does not affect the essential 
features of the model but leads to a considerable simpli- 
fication of the mathematical apparatus. 

The starting point is a many-body Hamiltonian with 
pairing residual interactions. This Hamiltonian depends 
on some time-dependent collective parameters q{t) - 
{qi{t)] (/ = l,...n), such as the internuclear distances 
between atoms or nuclei: 

H{t) = ^ £k[q{t)]{alak+a^aj,)-G ^ fl^fl^fl,•fl;.(l) 
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Here, are single-particle energies of the molecular 
potential, and denote operators for creating and 
destroying a particle in the state k, respectively. The 
state characterized by a bar signifies the time-reversed 
partner of a pair The pairing correlation arise from 
the short range interaction between fermions moving in 
time-reversed orbits. The essential feature of the pair- 
ing correction can be described in terms of a constant 
pairing interaction G acting between a given number of 
particles. In this paper, the sum over pairs generally runs 
over the index k. Because the pairing equations diverge 
for an infinite number of levels, a limited number of lev- 
els are used in the calculation, that is levels above and 
below the Fermi energy Ef. 

Using quasiparticle creation and annihilation opera- 
tors 



ak(y) - 



+ 

k(y) 



a 



Uk{y)ak - Vk(y)Cl-p. ; 

Uk(y)ai + Vk(y)ap, 

Uk(y)at - K(yfk\ 
Uk(y)al + Vk(y)''k', 



(2) 



it is possible to construct some interactions able to break 
a Cooper pair when the system traverses a avoided 
crossing region. The parameters Vk(y) and Uk{y) are occu- 
pation and vacancy amplitudes, respectively, for a pair 
occupying the single-particle level k of the configuration 
(y). The seniority-zero configuration is labeled with 
(y) - (0) and the seniority-two configuration by a pair 
of indexes denoting the levels blocked by the unpaired 
fermions (y) = {ij). The three situations plotted in Fig. 
[T]can be modeled within products of such creation and 
annihilation operators acting on Bogoliubov wave func- 
tions. In the plotflja), the Cooper pairs remain on the 
adiabatic levels e, and ej after the passage through the 
avoided crossing region, in Fig. [I](b) the pair destruc- 
tion is illustrated, while in Fig. ]l\{c) two fermions gen- 
erate a pair after the passage through an avoided cross- 
ing region. Formally, to describe these three situations, 
an interaction in the avoided crossing can be postulated 
as follows: 
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where hij is the interaction between the levels. The 
form of the perturbation ^ was postulated in Ref. ifllll 
and was successfully used to generalize the Landau- 
Zener effect in seniority-one superfluid systems. Act- 
ing on a suited Bogoliubov wave function, the product 




Figure 1 : Ideal avoided crossing regions between two adiabatic lev- 
els e, and ej. Thi'ee possible transitions between configurations in 
an avoided crossing region in the superfluid model are displayed: (a) 
The configuration remains unchanged after the passage through the 
avoided crossing region, (b) A parr is broken, (c) A pair is created. 



over k transforms the seniority-two configuration in the 
seniority-zero one in the case of the first term in the left 
hand of Eq. and vice-versa in the case of the second 
term. If the product ck,(0)Q'J(0) acts on a seniority-zero 
function, then it annihilates a pair and creates two un- 
paired fermions in states / and j. If the product a^QjO-t 
acts on a seniority-two wave function, then it creates 
a pair distributed on both orbitals / and j. In order to 
obtain the equations of motion, we shall start from the 
variational principle taking the following energy func- 
tional 
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and by assuming the many-body state formally ex- 
panded as a superposition of time dependent BCS 
seniority-zero and seniority-two adiabatic wave func- 
tions 
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where cq and cji are amplitudes of the two kind of con- 
figurations. Here, A is the chemical potential, and N is 
the particle number operator. To minimize this func- 
tional, the expression (|4| is derived with respect the 
independent variables Vk(0), Vk(ji), cq, cji, together with 
their complex conjugates, and the resulting equations 
are set to zero. Eventually, eight coupled-channel equa- 
tions are obtained: 

ifipkiO) = '<k{0)^Q - 1(1(0)^0', 

i^PkUD = Kk(ji)i^), - 4ui)^Ji'- 
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where the partial derivatives with respect the time are 
denoted by a dot. The sums are restricted by the condi- 
tions j + I, m + n, m + j, and n + I. Ey are exactly the 
expected values of the Hamiltonian ([T]) for the seniority- 
zero or seniority-two configurations: 
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and Tkiy) are energy terms associated to single-particle 
states: 
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The following notations are used in Eqs. 
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where pk(y) are single-particle densities and Kk{y) are 
pairing moment components. Py denote the probabil- 
ities to find the system in the configurations y. S yy are 
moment components between two configurations y and 
y and have the property | Syy PyPy. Ay is the gap 
parameter. The values of pk(y) and Py are reals. The par- 
ticle number conservation conditions 2 Y^k Pk(0) - 2A^, 
2 llki=j,iPk{ji) ^2N -2 and Po + Hjj^jPji = 1 ai'e ful- 
filled by Eqs. ©. 



3. Application to fission processes 

A direct application of the system ^ is related to the 
pair breaking and the odd-even structure in fission frag- 
ment yields. In fission, it is considered that the paired 
configuration is preserved until an interaction breaks 
some pairs in combination with the existence of a suffi- 
cient high excitation energy. The odd-even structure in 
fission is explained usually within statistical arguments, 
as for example in Refs. |13L . Alternatively, the prob- 
ability to break a pair can be determined dynamically by 
taking into account only the interaction available in the 
avoided crossing regions within the present model. 

To solve the pair breaking equations the varia- 
tions of single-particle energies ek together with per- 
turbations hjj must be supplied. The simplest way to 
obtain these evolutions is to consider a time-dependent 
mean field in which the nucleons move independently. 
In most usual treatments of nuclear fission, the whole 
nuclear system is characterized by some collective vari- 
ables, which determine approximately the behavior of 
many other intrinsic variables. The generalized coor- 
dinates vary in time leading to a split of the nuclear 
system. The basic ingredient in such an analysis is a 
nuclear shape parametrization. In the following treat- 
ment, a nuclear shape parametrization is given by two 
ellipsoids of different sizes smoothly joined by a third 
surface obtained iii|] by the rotation of a circle around 
the axis of symmetry. Five degrees of freedom charac- 
terize this parametrization: the elongation given by the 
inter-nuclear distance R between the centers of the el- 
lipsoids, the two deformations of the nascent fragments, 
the mass asymmetry and the necking parameter. Due to 
the axial symmetry, the good quantum numbers are the 
projections of the intrinsic spin Q.. 

As specified in Ref. [15], first of all, a calculation of 
the fission trajectory in our five-dimensional configura- 
tion space, beginning with the ground-state of the sys- 
tem up to the exit point of the barrier must be performed. 
This can be done by minimizing the action integral. For 
this purpose, two ingredients are required: the defor- 
mation energy V and the tensor of the effective mass. 
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Figure 2: (a) Proton level scheme as function of the inter-nuclear dis- 
tance R between the nascent fragments along the minimal action path 
for the ^^"'U fission. Only N=60 single-particle energies above and 
below the Fermi level are plotted. The Fermi level Ef is displayed 
with a thick full curve. The selected £1=1/2 levels (ei, ...en) are plot- 
ted with thick dot-dashed lines. Avoided level crossing regions are 
marked with circles. The maximum number of major shells used in 
calculations is 12. (b) The fission barrier as function of the inter- 
nuclear distance R. Nuclear shapes corresponding to several special 
configurations are inserted in the plot. The values of the intemuclear 
distances are also marked. At R » 4 fm, the ground-state configura- 
tion of the parent is found. 



The deformation energy was obtained in the frame of 
the microscopic-macroscopic method 1 16] by summing 
the liquid drop energy with the shell and the pairing 
corrections. The macroscopic energy is obtained in the 
framework of the Yukawa plus exponential model IitIi 
extended for binary systems with different charge densi- 
ties [18]. The Strutinsky microscopic corrections were 
computed on the basis of the Woods-Saxon superasym- 
metric two center shell model. The effective mass is 
computed within the cranking approximation |19|] . Af- 
ter minimization, the dependences between the general- 
ized coordinates qt (/ = 1, ...5) in the region comprised 
between the parent ground state configuration and the 
exit point of the external fission barrier supply the least 
action trajectory. The ground-state corresponds to the 
lowest deformation energy in the first well. The least 
action trajectory is obtained within a numerical method. 
Details about the numerical procedure of minimization 
and about the model can be found in Refs. [10, 1 1] and 
references therein. The resulting ^^*U fission barrier is 
plotted on Fig. |2](b) as function of the distance between 
the centers of the nascent fragments R. Some nuclear 
shapes obtained along the minimal action trajectory are 
inserted in the plot. A realistic proton level scheme 
along the least action trajectory is obtained within the 
superasymmetric Woods-Saxon two-center shell model 
[11]. This model gives the single particle level diagrams 
by diagonalizing a Woods-Saxon potential, corrected 
within spin-orbit and Coulomb terms, in the analytic 
eigenvalue basis of the two center semi-symmetric har- 
monic model [20, 21]. Other recipes to obtain the level 
scheme are related to the molecular orbital method ||22|] . 
The proton level scheme is displayed in Fig. |2](a). 

The Landau-Zener effect is produced only in the 
workspace spanned by levels characterized by the same 
good quantum numbers. As mentionned, due to the ax- 
ial symmetry, the good quantum numbers are the pro- 
jections of the intrinsic spin Q. Therefore, among the 

states in the region near the Fermi surface, the 17 
levels with spin projection Q=l/2 are selected. These 
levels are plotted with thick dot-dashed lines in Fig. |2] 
(a). Because a pair creation or annihilation is consid- 
ered to be possible only between adjacent levels, 16 
seniority-two configurations are constructed. The ini- 
tial values of quantitiesdSJ are obtained by solving the 
BCS equations for the ground states of all y configura- 
tions. The average excitation energy of the seniority- 
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where is the value 



zero state is computed as in Ref. 
relation E* = (£0 - E°) + E* 
of the lowest energy state of any deformation calcu- 
lated within the BCS approach and £() is obtained with 
Eqs. ©. Here, £* is the dissipation obtained within 
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Figure 3: (a) The average excitation energy E* as function of the elon- 
gation R. The inter-nuclear velocities R are 10", 10^ and 10*' m/s for 
the dashed, full and dot-dashed lines, respectively, (b) The probabil- 
ity Podd of a seniority-two state with respect to the elongation R. The 
same line types and inter-nuclear velocities as in panel (a) are used, 
(c) The probabilities to obtain a seniority-zero state Pq and a seniority- 
two state Pgjj as function of the excitation energy E' of the fragments 
at the elongation R=2Q fm. 



the same formalism for the even neutron subsystem. 
The probabihty to obtain a seniority-two state is sim- 
ply Podd = 1 - ^'o- The values of E* and Podd are 
plotted in Fig. |3](a) and (b) as function of the elonga- 
tion for some inter-nuclear velocities R. In connection 
with the shape of the barrier displayed in Fig. |2l it can 
be deduced that the larger part of the odd-even yield is 
formed during the penetration of the second barrier and 
the excitation energy increases merely in the same re- 
gion. Different constant values of the inter-nuclear ve- 
locity R ranging from 10"* to 10*' m/s were tested. These 
values can be translated in a time to penetrate the bar- 
rier ranging in the interval [1.4 x 10"'**, 1.4 x 10"^"] s. 
In Fig. [3](c), the dependences of Pq and P„dd versus E* 
are displayed in the selected velocity domain. The re- 
sults exhibit a clear decrease of Podd as function of E*. 
It is interesting to note that at zero excitation energy, the 
probability to find the system in a seniority-two state is 
practically one. In cold fission, at very low excitation 
energies of the fragments, the odd-even yields are al- 
ways larger than the even-even ones ll25[ [26 
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even-even fragmentation dominates at larger excitation 
energies of the fragments, above 4-6 MeV. It is a very 
strange phenomenon because in cold processes the sys- 
tem doesn't possess enough energy to break a pair and 
because the penetrability is hindered for odd-systems 
due to the specialization energies associated to the two 
unpaired nucleons. Up to now, the statistical explana- 
tion of this phenomenon involved some modifications 
of the level densities for odd-even and even-even par- 
titions lfl4ll by according them within the deformations 



of the fragments as function of the excitation energies 
and the shell effects. However, if one assumes that the 
odd-even effect in the fission fragments distribution is 
strongly correlated to the seniority-two state probability, 
this phenomenon can be alternatively explained by solv- 
ing the coupled-channel system of time-dependent pair 
breaking equations as evidenced above. In this work, 
only the Q=l/2 subspace of the proton level diagram is 
treated, but the same formalism can be applied to other 
subspaces. 

In conclusion, a new set of time-dependent coupled 
channel equations derived from the variational principle 
is proposed to determine dynamically the mixing be- 
tween seniority-zero and seniority- two configurations. 
The essential idea is that the configuration mixing is 
managed under the action of some inherent low ly- 
ing time dependent excitations produced in the avoided 
crossing regions. These equations were used to explain 
the odd-even effect in cold fission processes. Only the 
radial coupling was used in the analysis, but it is pos- 
sible to extend the equations to take into account even 
the Coriolis coupling, as in Ref. [28]. The main trends 
concerning the dependence of the odd-even effect in 
fragments yields versus the fragments excitation energy 
were reproduced. The formalism can be adjusted for 
other types of processes. In this respect, the value of the 
interaction hyy can also be the magnitude of other kind 
of interactions between different configurations and the 
collective variables could be the amplitude of an ex- 
ternal applied field as encountered in the field of con- 
densed matter The essential idea is the inclusion in the 
energy functional of a time dependent low perturbation 
by the mean of quasiparticle operators. 
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